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SECTION  1 


INTRODUCTION 

The  purpose  of  this  investigation  is  to  determine  the  perfor¬ 
mance  payoff  of  post  stall  maneuvering  (PSM>,  sometimes  called 
supermaneuverability ,  and  thrust  vectoring.  FSM  involves 
maneuvering  at  high  angles  of  attack,  angles  greater  than  that  for 
max 'mum  lift  coefficient.  There  are  two  potential  PSM  advantages. 
The  first  is  the  ability  to  rapidly  change  the  flight  path  of  the 
vehicle.  The  second  is  to  obtain  the  larger  pointing  angles  which 
is  important  from  a  weapon's  standpoint.  Pointing  will  not  be 
addressed  in  this  effort,  but  clearly  it  has  an  advantage  at  any 
speed . 

In  the  past,  conventional  maneuvers  were  restricted  to  less 
than  maximum  lift  angles  of  attack.  If  PSM  has  any  advantage,  it 
must  be  able  to  develop  velocity  rotational  rates  greater  than 
that  for  conventional  maneuvers.  One  way  of  addressing  the  PSM 
issue  is  by  solving  instantaneous  turning  performance  and  minimum 
time  to  turn  problems.  That  is  the  approach  that  will  be  taken 
here.  The  analysis  will  be  limited  to  subsonic  speeds.  The 
reason  for  this  is  that  past  performance  studies  have  shown  that 
if  there  is  a  performance  payoff  for  PSM,  then  it  occurs  at  low 
speeds;  that  is  speeds  below  the  corner  velocity. 

The  first  group  of  problems  will  be  instantaneous  turning 
performance  problems.  Thrust  vectoring  will  be  included.  The 
solution  to  this  problem  yields  insight  into  both  the  effect  of 
thrust  vectoring  and  high  angle  of  attack  maneuvers  along  the 
trajectory.  The  last  problem  is  the  variable  speed  optimal  con¬ 
trol  problem  for  minimum  time  turns.  The  significance  of  this 
problem  is  that  the  effect  of  high  angle  of  attack  on  drag  as  well 
as  lift  during  accelerated  maneuvers  can  be  addressed. 


An  important  difficulty  with  high  angle  of  attack  maneuvers 
he  ability  to  control  the  vehicle.  This  issue  will  not  be 
essed  in  this  study.  It  will  be  assumed  that  the  vehicle  can 
ontrolled  anywhere  between  zero  and  ninety  degrees  angle  of 
ck.  This  assumption  should  be  addressed  on  its  own  since  it 
tes  to  a  critical  problem  area  in  achieving  PSM  capability. 

Several  papers  (References  1-5)  have  focused  attention  on 
mum  time  turns  for  high  speed  aircraft.  These  efforts  include 
two  and  three  dimensional  trajectories.  The  former  cor- 
ond  to  constant  altitude  or  vertical  plane  maneuvers,  either 
t-S  or  half  loops.  In  the  next  section,  the  equations  of 
on  and  the  control  constraints  are  presented. 


2 


tu  -  —  L 
max  v  max 


(50) 


Lmax  *  «SCL  'W 

max 


(51) 


The  global  maximum  of  the  turning  rate  is  obtained  by  maxi- 
ing  simultaneously  with  respect  to  both  the  angle  of  attack  and 
thrust  vector  angle.  Recall  Equation  (15) 


ui  *  —  CT  sin(<x+S  )+L  3 
V  v 


( 15) 


optimum  is  obtained  by  maximizing  both  terms.  Thus,  the  angle 

attack  is  that  for  Cr  and  the  thrust  vector  angle  is  the 

max 

plement  of  the  angle  of  attack.  Hence 


vjj  =  -  (T  +L  ) 
max  „  max 

v 


(52) 


s  result  can  also  be  determined  rigorously.  In  order  for  to 

aaximized,  it  is  necessary  that  the  following  equations  be 
Lsf ied 


sr  *  t’co*<«+sv>«-  §£- 


(53) 


"tt —  =  T  cos(a+S  >  =  0 
dSv  v 


(54) 


n  the  latter,  the  sum  is  ninety  degrees.  Substitution  into  the 
5 1  equation  shows  that 
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The  optimum  angle  of  attack  for  a  given  Sv  is  obtained  from 
Equations  (15)  and  (24) 

Tcos(a+Sv>+qSCL  =  Tcos(a+Sv)+qS(C^  cosa-C^sina)  =  0  (45) 

a  a 


The  optimum  angle  of  attack  is  greater  the  °c  - Ldl i *  optimum 

k 

angle  of  attack,  a  ,  is  the  solution  of 


A  Tcos&v+qSCN 

tana  =  _  a 

TsinS  +qSC^ 


(46) 


This  equation  agrees  with  Equation  (26)  when  Sv  =  0.  At  low 
speeds  where  q  is  small. 


tana  ■  cots 


The  solution  is 


*  _  1  c 

a  "  2  Sv 


(48) 


Since  Tsin(a  +Sv>  dominates,  from  Equations  (14)  and  (15)  maximum 
hi  is  approximately 


g 

ut  *  —  T 

max  v 


(49) 


At  high  speeds  where  q  is  large.  Equation  (45)  shows  that  CL 


must  decrease  and  a  therefore  approaches  astail*  ^  dominates 


U  is  achieved  at  the  same  »gtall*  ^as  been  observed  in  ex- 
erimental  aerodynamics  that  varies.  At  low  Mach  numbers 

here  the  thrust  dominates,  turning  rate  is  approximately 


9 

w  -  -  T  <44) 

V 

nd  is  insensitive  to  aerodynamic  variations.  At  high  Mach  num- 
ers  where  the  lift  dominates,  the  turning  rate  varies  in  a  linear 
anner  with  respect  to  the  maximum  value  of  the  normal  force 
oef f icient . 


The  explanation  for  this  is  as  follows.  At  low  speeds,  a  is 
arge  and  approximately  ninety  degrees.  Equation  (28)  shows  that 
he  maximum  turning  rate  is  insensitive  to  variations  in  C^.  At 

igh  speeds.  Equation  (43)  shows  that  the  maximum  turning  rate  is 
roportional  to  the  normal  force. 

The  optimum  angle  of  attack  is  presented  in  Figure  13.  At 

aw  speeds.  Equation  i26)  shows  that  a*  decreases  with  increasing 

,.  .  At  high  speeds,  a*  approaches  which  we  assumed  to 

max 

s  independent  of  maximum  for  this  steady  state  aerodynamic 


irust  Angle  Variations 

The  analysis  of  the  influence  of  Sv  is  next.  In  Figures  14 

id  15  are  presented  the  maximum  turning  rate  and  the  correspond- 
ig  angle  of  attack,  respectively.  The  thrust  angles  vary  from 
*ro  to  sixty  degrees  in  ten  degree  increments.  The  maximum 
irning  rate  is  relatively  insensitive  to  the  thrust  vector  angle. 


14 


where 


N 


qscN 

_ A 

w 


(42) 


,  ,  /  2  .2,2 

At  high  speeds,  CL  =  0,  N  *  N  tana,N  >>T  N  »T  N  ,  hence 

a  a  a 


w 


max 


g 

*  —  N  cosoc 
V 


A 


(43) 


The  maximum  turning  rate  is  relatively  insensitive  to  thrust 
variations. 


The  sensitivity  of  the  optimum  angle  of  attack  to  variations 

in  the  thrust  is  presented  in  Figure  11.  At  low  speeds  where  a* 
is  large.  Equation  (26)  shows  that  tana  varies  linearly  with  the 
thrust.  At  high  speeds.  Equation  (23)  shows  that  approaches 


zero. 


Hence  a 


A 


approaches  «#tall 


which  is  independent  of  the 


thrust . 


Aerodynamic  Variations,  Sv  *  0 


A  sensitivity  analysis  of  the  influence  of  the  maximum  value 

of  normal  coefficient  on  the  maximum  turning  rate  was  conducted. 

The  variation  in  the  aerodynamic  normal  force  coefficient  is  a 

constant  multiple  of  the  values  defined  earlier.  C„  data  are 

wnom 

the  nominal  values  used  earlier.  Maximum  turning  rate  is 

presented  in  Figure  12.  The  value  of  C^/C^  =  2.0  represents  an 

nom 

extreme  change  in  steady  state  aerodynamics.  In  unsteady  flow,  it 
may  be  possible  to  exceed  this  value.  It  is  assumed  that  maximum 
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(38) 


V*  a  -JT _ 

CRIT  pSCL 


The  critical  speed  is  independent  of  a.  For  high  angles  of  at¬ 
tack,  sina  increases,  decreases  and  Equation  (37)  shows  the 
critical  speed  increases.  Figure  9  shows  that  the  critical  speed 
is  approximately  0.15  Mach  number.  The  equilibrium  speed  was 
approximately  0.2  Mach  number.  The  minimum  turning  rate  is  ob¬ 
tained  by  substituting  Equation  (36)  into  Equation  (15) 


2g 

—  T  sina 
V 


(39) 


Thrust  Magnitude  Variations,  Sv  *  0 

A  sensitivity  analysis  with  respect  to  maximum  thrust  varia¬ 
tions  was  conducted.  For  an  altitude  of  20,000  feet,  maximum 
turning  rate  is  presented  in  Figure  10  for  different  thrust 

levels.  T _ _  corresponds  to  the  nominal  values  used  in  the  pre- 

vlous  calculations.  At  low  speeds  the  maximum  turning  rate  is 
nearly  linear  with  respect  to  the  thrust.  At  high  speeds  the 
turning  rate  is  essentially  insensitive  to  the  thrust  since  the 
lift  dominates. 


This  behavior  can  be  explained  by  examining  Equation  (28) 
At  low  speeds,  aA  is  large,  thus  approximately 


g  - 

-  T 
V 


(40) 


Maximum  turning  rate  varies  linearly  with  the  thrust.  The  situa¬ 
tion  at  high  speeds  is  explained  as  follows.  Combining  Equation 
(15),  (21),  and  (26)  gives  after  some  manipulation 


u» 


max 


gcosa 

- —  <T 

VN 


2  ,2 

+N  +T  N  ) 
a 


(41) 
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thrust.  This  is  not  to  be  confused  with  absolute  minimum  turning 
rate  which  corresponds  to  zero  lift  and  zero  thrust. 

For  a  given  angle  of  attack.  Equation  (15)  and  Figure  7  shows 
that  u  varies  with  respect  to  V  in  the  following  nanner 

A 

u»  =  -  +  BV  (33) 

V 


where 


,  gpSCr 

A  =  gT  sina,  B  *  -  (34) 

2H 


A  and  B  are  invariant  with  respect  to  the  speed  if  the  thrust  does 
not  vary  with  respect  to  the  speed.  If  the  thrust  can  be  treated 
as  a  constant,  then  the  minimum  turning  rate  is  determined  from 
the  minimum  of  Equation  (33) 


A 

-  —  +  B  *  0  (35) 

V2 


or 


Tsina  =  L 

The  critical  speed,  VCRIT,  is  the  solution  of 

.. z  2Tslna 

v  CRIT  pSCL 


(36) 


(37) 


For  low  thrust,  VCRIT 
approximately  a  and 


is  low.  For  small  angles  of  attack,  sina  is 
is  approximately  a.  Thus  Equation  (37) 


is  approximately 
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(32) 


Thus  the  angle  of  attack  for  maximum  turning  rate  and  V  =  0, 

y  =  0  are  the  same  whenever  the  maximum  occurs  at  an  angle  of 

attack  where  CN  is  constant.  Refer  to  the  intersection  between 

the  two  angle  of  attack  curves  as  the  equilibrium  speed.  The 

reason  for  calling  it  this  is  that  the  speed  will  not  change  once 

this  point  is  reached  because  \?  =  0.  The  only  way  the  speed  can 

decrease  is  if  the  angle  of  attack  is  increased  creating  more  drag 

or  the  thrust  is  reduced.  It  can  be  concluded  that  if  the  initial 

speed  is  greater  than  the  equilibrium  speed  and  C^  is  equal  to 

C„  the  the  vehicle  declerates  to  the  equilibrium  speed  and  then 
Nmax 

remains  there.  If  the  flight  path  angle  is  not  zeror  then  from 
Equations  (1)  and  (4)  the  equilibrium  speed  is  less  (greater)  than 
the  zero  flight  path  case  if  the  flight  path  angle  is  positive 
(negative) . 

At  higher  speeds  the  structural  limit  comes  into  play.  The 
speed  at  which  this  occurs  increases  with  increasing  altitude 
since  the  dynamic  pressure  decreases.  Only  the  structural  limit 
for  ten  thousand  feet  is  shown  since  the  corner  speed  at  twenty 
thousand  feet  is  greater  than  a  Mach  number  of  0.6.  The  angle  of 
attack  for  V  *  0,  y  =  0  is  less  than  that  for  maximum  turning  rate 
since  the  drag  in  Equation  (4)  increases  with  increasing  speed. 

In  order  to  satisfy  V  =  0,  the  angle  of  attack  must  decrease  with 
increasing  speed.  In  Figure  9  are  presented  turning  rates  cor¬ 
responding  to  the  maximum,  V  -  0  for  y  =  0,  and  the  structural 
limit  at  ten  and  twenty  thousand  feet.  As  the  altitude  Increases, 
the  thrust  decreases  for  a  given  Mach  number  and  the  dynamic 
pressure  decreases  al3o.  Consequently,  optimum  u>  decreases. 

Critical  Speed,  6v  *  0 

The  critical  speed  is  that  speed  where  the  turning  rate  is 
minimum  with  respect  to  speed  for  a  given  angle  of  attack  and 
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At  low  speeds,  q  is  small  and  a  approaches  ninety  degrees.  At 
high  speeds,  q  is  large  and  <x*  approaches  agtan*  conclusion 

is  deduced  from  comparison  with  Equation  (24).  Equation  (23)  can 
be  rewritten  as  follows 


Tcosa+L  = —  *  0 


(27) 


Solving  for  L  and  substituting  into  Equation  (15)  gives 


9  *  C! 

-  T'  (  sina  -  — * 

V  C, 


A 

cosa  ) 


(28) 


In  Figure  8  is  illustrated  the  optimum  angle  of  attack  for 
altitudes  of  ten  and  twenty  thousand  feet.  Also  included,  are  the 
angles  of  attack  corresponding  to  zero  longitudinal  acceleration 
and  a  load  factor  limit  equal  to  7.33.  At  speeds  below  the  corner 
speed,  the  structural  limit  can  not  be  reached.  Also,  at 
low  speeds  the  maximum  turning  rate  and  the  turning  rate  at  V  =  0 
and  zero  flight  path  angle  are  the  same.  The  reason  for  this  is 
the  angle  of  attack  is  higher  than  that  value  where  reaches  its 
maximum  value.  Since  maximum  turning  rate  corresponds  to  maximum 
,  from  Equations  (5)  and  (21) 


Fn  *  T  sina+N  cosa 


(29) 


'  >  a»j, 

rr~~  *  T  cosa-N  sina+  —  cosa 
da  da 


(30) 


But  the  last  term  is  zero  since  is  constant.  Comparison  with 
Equation  (4) 


T  cosa-N  sina 


(31) 


shows  that 


location  ot  the  maximum  value  are  sensitive  to  the  Mach  number. 

In  Figure  7,  u  is  presented  as  a  function  ot  Mach  number  for 
constant  values  of  the  angle  of  attack.  The  latter  varies  from  0 
to  90  degrees  in  5  degree  increments.  For  constant  a,  u>  decreases 
with  increasing  Mach  number  until  the  minimum  is  reached.  This 
minimum  will  be  referred  to  as  the  critical  speed.  From  there  on 
u>  increases  through  the  subsonic  region. 


For  a  given  Mach  number,  there  is  an  angle  of  attack  that 
gives  maximum  u>.  At  low  Mach  numbers  the  angle  of  attack  is  high 
and  goes  to  90  degrees  as  the  Mach  number  goes  to  zero.  At  high 
Mach  numbers  the  optimum  angle  of  attack  approaches  that  for 
aSTALL*  •^his  behavior  can  be  seen  by  inspection  of  Equations  (14) 
and  (15).  Optimization  of  Equation  (14)  with  respect  to  the  angle 
of  attack  requires  that 


3L 


Tcoscx+  —  =  Tcosa+qSC^ 


(23) 


Since  the  first  term  is  positive  the  second  term  must  be  negative. 
This  can  occur  only  for  angles  of  attack  greater  than  that  for 
aSTALL '  From  Equation  (21) 


cosa-C^sina 

(X 


(24) 


Equation  (23)  becomes 


Tcosa+qS(CN  cosa-C^sina)  =0  (25) 

(X 

The  optimum  angle  of  attack,  a*,  for  a  given  speed  and  altitude  is 
determined  from 


tana 


A 


T+qSC 

_ jx 

qSCN  ' 


(26) 
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Cjj  *  CLcosa+CDsina 


(17) 


=  Cjjcosa-C^sina  (18) 

The  steady  state  data  correspond  to  subsonic  Mach  numbers.  It 
will  be  assumed  throughout  this  effort  that  the  Mach  numbers  are 
sufficiently  low  that  Mach  number  effects  can  be  neglected.  The 
coefficients  that  correspond  to  .the  data  of  Figures  2  and  3  are 
presented  in  Figure  5.  From  these  data,  approximate  relations  can 
be  deduced  for  use  in  performance  prediction 

CA  *  0  (19) 

*  0.16+3.44<x  0<  a  <0.44  rad 

«  -1. 72+12. 04a-9. 86a2  0.44<  a  <0.61  rad  (20) 

=  1.96  0.61<  a 

The  relations  for  the  lift  and  drag  coefficients  reduce  to 


CL 

=  C^cosa 

(21) 

CD 

=  Casino, 

(22) 

These  equations  are  exact  if  is  exactly  zero.  They  are  in  good 
agreement  with  the  experimental  data  as  can  be  seen  by  examination 
of  Figures  2  and  3.  The  approximate  relations  are  used  in  all 
subsequent  calculations. 

Angle  of  Attack  Variations,  Sv  =  0 


For  no  thrust  vectoring,  $v  =  0.  For  a  given  Mach  number, 
T/W,  and  altitude,  the  turning  rate  looks  like  that  illustrated  in 
Figure  6.  The  turning  rate  has  a  maximum  somewhere  between  0  and 
90  degrees  angle  of  attack.  Both  the  shape  of  the  curve  and  the 
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The  maximum  of  y  corresponds  to  a  =  0  and  maximum  w.  Minimum  y 
corresponds  to  a  *  it  and  maximum  u.  Thus  in  both  cases,  the 
limits  result  from  maximum  w.  From  Equation  (15)  it  follows  that 
maximum  u>  results  from  maximizing  F^.  Clearly  maximum  thrust 
maximizes  F^.  The  rest  of  the  section  will  concentrate  on  maxi¬ 
mizing  the  turning  rate  function  w. 

Aerodynamic  Relations 

Typical  aerodynamic  data  from  Reference  6  are  illustrated  in 
Figures  2  and  3,  and  correspond  to  the  F15  which  is  a  high  perfor¬ 
mance  fighter.  The  data  represent  steady  state  aerodynamics.  Non 
steady,  or  time  dependent,  aerodynamic  lift  data  are  superimposed 
on  top  of  the  steady  state  data  in  Figure  4.  The  non  steady  data 
are  only  an  illustration  of  a  phenomena  called  dynamic  overshoot 
and  undershoot.  When  the  rate  of  change  of  the  angle  of  attack  is 
very  large,  the  available  lift  first  exceeds  the  steady  state 
lift.  After  the  stall  angle  of  attack  is  reached,  there  is  a 
rapid  transition  from  overshoot  to  undershot  where  the  lift 
decreases  with  angle  of  attack.  When  the  angle  of  attack  is 
rapidly  reduced,  the  dynamic  lift  drops  below  the  steady  state 
values.  The  overshoot  and  undershoot  aerodynamics  are  not  com¬ 
pletely  understood  to  say  the  least,  therefore,  there  is 
uncertainty  in  their  prediction.  There  is  some  understanding  of 
the  variables  involved,  however.  The  overshoot  is  primarily  a 
function  of  a  and  the  final  value  of  a.  The  transition  from 
overshoot  to  undershoot  is  time  dependent  and  the  lift  may  be 
asymmetric,  one  wing  may  stall  before  the  other.  The  undershoot 
depends  on  the  value  of  maximum  in  overshoot  and  is  also  time 
dependent.  The  drag  also  experiences  variations  from  the  steady 
state  values.  The  drag  may  lag  the  non  steady  lift.  For  the  rest 
of  this  effort,  we  will  consider  only  steady  state  aerodynamics . 

The  normal  and  axial  force  coefficients  are  related  to  the 
lift  and  drag  coefficients  and  the  angle  of  attack  in  the  follow¬ 
ing  way 


IU  . 
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SECTION  III 


INSTANTANEOUS  TURNING  PERFORMANCE 

Instantaneous  turning  performance  corresponds  to  point  per¬ 
formance  analysis.  Each  point  in  the  flight  envelope  can  be 
studied  independently  of  every  other  point.  In  Section  IV  where 
continuous  turning  performance  is  studied,  it  is  proven  that  for 
minimum  time  turns  the  optimal  trajectory  is  in  the  vertical 
plane.  Therefore,  this  section  will  consider  only  vertical  plane 
trajectories.  The  appropriate  turning  rate  equation  is  the  time 
rate  of  change  of  the  flight  path  angle 

g 

y  -  -<Fmcoso-cosy)  (13) 

V  w 


where 


/  / 

FN  *  T  sin(a+6v)+L  (14) 

O 

The  bank  angle  o  is  zero  or  180  degrees  (0  corresponds  to  a  pull 

O 

up  and  180  corresponds  to  a  pull  down).  Introduce  the  turning 
rate  function  w  defined  by 


g 


The  actual  flight  path  angular  rate  is  related  to  u  as  follows 

g 

Y  =  wcoso-  —  cosy  (16) 

V 

At  any  point  on  the  trajectory,  that  is  given  the  speed,  flight 
path  angle,  and  the  altitude,  the  objective  is  to  determine  the 
propulsive  and  aerodynamic  controls  that  optimize  the  flight  path 
angular  rate.  This  includes  both  maximizing  and  minimizing  y. 
This  is  accomplished  by  optimizing  the  turning  rate  function  u>. 
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Z  =  VsiJTf 


(9) 


The  change  in  the  mass  will  be  neglected.  The  constraints  for 
these  equations  are 


n  <>  n 


(10) 


T  <  T  <  T 
min  -  1  -  max 


(11) 


For  this  study  minimum  thrust  is  zero.  There  are  four  controls 
for  this  situation.  The  propulsion  controls  are  the  thrust  mag¬ 
nitude.  T,  and  thrust  vector  angle,  V  The  aerodynamic  controls 
are  the  aerodynamic  angle  of  attack,  a,  and  bank  angle,  o.  Since 
the  study  has  application  to  post  stall  maneuvering,  the  angle  of 
attack  is  limited  only  by  the  maximum  value  of  the  structural  load 
factor  and  ninety  degrees.  In  previous  studies  it  was  limited  to 
the  stall  angle  of  attack,  If  the  speed  is  greater  than 

the  corner  speed,  then  the  angle  of  attack  is  limited  by  the 
maximum  value  of  the  load  factor.  The  corner  speed  occurs  at  the 

juncture  between  the  Cr  and  n.„„  limits.  It  is  determined  from 

u  max 

max 


7  pSCL  VC  *  "max" 

max 


(12) 


In  the  next  section,  instantaneous  turning  performance 
problems  are  addressed. 


SECTION  II 


PROBLEM  FORMULATION 

The  coordinate  system  selected  for  the  trajectory  problems  is 
a  velocity  coordinate  frame.  The  definition  of  the  velocity 
vector  using  a  velocity  coordinate  frame  includes  the  magnitude  of 
the  velocity,  the  heading  angle,  and  the  flight  path  angle.  The 
rate  equation  for  these  variables  are 


V 

-  g(FT-sinv) 

(1) 

• 

g 

Y 

*  — (FMcosa-cosy> 

V  N 

(2) 

• 

g 

<l> 

?  r  sina 
v  rN  siny 

(3) 

where 


ft 

=  T  cos(a+S  )-D 

V 

(4) 

fn 

*  T  sin(a+S  )+L 

V 

(5) 

f 

T 

=  T/W,  D  =  D/W,  L  *  L/W  =  n 

(6) 

The  aerodynamic  and  propulsive  forces  are  illustrated  in  Figure  1. 
The  thrust  is  a  function  of  throttle  setting,  Mach  number,  and 
altitude.  The  lift  and  drag  are  functions  of  Mach  number,  al¬ 
titude,  and  the  angle  of  attack.  The  position  of  the  vehicle  is 
determined  from  the  kinematic  relations 

x  *  Vcos-ycostJ)  (7) 


y  *  Vcosvsinto 


(8) 


which  occurs  at  astaii • 

In  Figures  16  and  17,  sustained  turning  rate  and  the  cor¬ 
responding  angle  of  attack  are  presented.  Comparison  of  Figures 
14  and  16  shows  that  the  angle  of  attack  at  low  speeds  is  the  same 
for  maximal  turning  ra^e.  Since  is  nearly  constant ,  Equation 

(45)  reduces  to  Equation  (4).  Thus  for  zero  flight  path  angle, 
the  solution  for  maximum  and  sustained  turns  are  the  same.  At 
high  speeds,  maximum  turning  occurs  simultaneously  with  lon¬ 
gitudinal  deceleration.  Thus,  thrust  vectoring  will  result  in 
faster  decelerations  relative  to  the  zero  vector  angle  case. 

There  is  no  equilibrium  speed  if  thrust  vectoring  is  used  because 
for  maximum  turning  rate  a+6v  is  ninety  degrees  and  therefore  F^ 

is  the  negative  of  the  drag  and  $<0. 

P3M  Versus  Thrust  Vectoring 

In  Figure  18,  maximum  turning  rate  is  presented  for  optimal 
PSM  and  thrust  vector  maneuvers.  In  general,  thrust  vectoring 
gives  slightly  higher  turning  rates.  The  corresponding  optimal 
angles  of  attack  are  presented  in  Figure  19.  For  thrust  vector¬ 
ing,  astan  i®  approximately  31.6  degrees  and  Sv  is  58.4  degrees. 

Longitudinal  decelerations  at  zero  flight  path  angle  are  presented 
in  Figure  20.  Thrust  vectoring  produces  greater  decelerations 
since  the  thrust  component  in  Equation  (4)  is  zero.  Figures  18 
and  20  illustrate  that  comparable  turing  performance  can  be 
achieved  through  PSM  maneuvers  or  thrust  vectoring. 


3ma&Lz 

Based  strictly  upon  steady  state  aerodynamics,  subsonic  Mach 
nuibers  less  than  0.6,  and  Instantaneous  point  performance  con¬ 
siderations,  very  high  angles  of  attack  (beyond  stall)  turns 
should  be  employed  only  at  low  speeds  Where  the  thrust  dominates 
the  lift.  The  situation  is  reversed  at  high  Mach  numbers  where 
the  angle  of  attack  approaches  the  stall  value.  Increasing  thrust 
gives  higher  turning  rates  particularly  at  low  Mach  numbers  and 
increasing  normal  force  coefficients  increases  the  turning  rates 
at  high  Mach  numbers.  The  equilibrium  Mach  number  (Mach  number 
where  maximum  turning  and  zero  rate  deceleration  at  zero  flight 
path  angle  are  equal)  is  approximately  0.2.  For  pull  ups  the 
equilibrium  speed  decreases  and  it  increases  for  pull  down 
maneuvers.  With  thrust  vectoring,  the  optimal  angle  of  attack  is 
astall  an<*  optimal  thrust  vector  angle  is  the  complement  of 

the  angle  of  attack.  This  holds  for  all  speeds.  Hence  post  stall 
maneuvering  (PSM)  is  not  necessary  for  turning  if  thrust  vectoring 
is  available.  Either  optimal  PSM  maneuvers  or  thrust  vectoring 
yield  comparable  turning  performance. 

Future  studies  should  include  the  following  aspects: 

1.  Unsteady  aerodynamics 

2.  Analysis  of  aircraft  pointing 

3.  Compare  turning  versus  pointing 

4.  Air  combat  analysis  MvN 
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SECTION  IV 


OPTIMAL  VARIABLE  SPEED  MANEUVERING  PERFORMANCE 


It  was  demonstrated  in  the  previous  section  on  instantaneous 
turning  performance  that  high  angle  of  attack  maneuvers  have  a 
performance  payoff  only  at  very  low  speeds.  This  occurs  because 
the  thrust  force  dominates  the  lift  at  low  speed  because  the 
dynamic  pressure  is  low.  As  the  speed  increases,  the  aerodynamic 
lift  dominates  over  the  thrust  and  the  maximum  turning  rates  tend 
to  be  near  astan*  this  section,  the  impact  of  both  aero¬ 
dynamic  and  thrust  control  will  be  addressed.  The  latter  Includes 
both  thrust  magnitude  and  vector  angle  variations.  It  will  be 
assumed  throughout  that  the  thrust  magnitude  does  not  vary  with 
angle  of  attack  or  thrust  angle.  The  equations  of  motion  are 
Equations  (1)  through  (3) 


*  gtFy-sinv) 

9 

3  -vFlIcoso-cos'y) 

V  * 

g  sino 

=  —  F  - 

y  N  cosy 

where  FT  and  F^  are  defined  by  Equations  (4)  through  (6) 

# 

FT  «  T  cos(a+$v) -D 
Fu  *  T  sin(a+6.)+L 

N  V 

T  «  T/W,  D  =  D/W,  L  *  L/W 


V 


y 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


The  altitude  variation  will  be  assumed  to  be  negligible. 

The  optimal  control  problem  is  the  minimum  time  to  turn  from 
V(o)  «  VQ,  y(o)  =  0  (56) 


to 


Yitf)  =  0  (57) 

where  VQ  is  specified.  If  the  motion  is  in  the  vertical  plane, 

then  the  final  flight  path  angle  can  be  taken  to  be  ±  «.  The 
final  speed  may  be  free  or  it  may  be  specified.  If  VQ  is  free, 

the  optimal  trajectory  is  free  to  go  to  low  speeds  if  this  is 
optimum.  The  final  heading  is  free. 

There  are  four  controls;  the  thrust  magnitude,  the  thrust 
vector  angle,  the  aerodynamic  angle  of  attack,  and  the  bank  angle. 
Optimal  controls  will  be  determined  for  all  four.  Then  the  case 
where  the  thrust  angle  is  fixed  and  set  equai  to  zero  will  be 
addressed.  The  contraints  on  the  structural  load  factor  and 
thrust  magnitude.  Relations  (10)  and  (11),  must  hold  everywhere 
along  the  trajectory.  The  problem  with  aerodynamic  control  at 
high  angles  of  attack  will  not  be  addressed.  It  will  be  assumed 
that  there  is  sufficient  control  to  fly  anywhere  between  zero  and 
ninety  degrees. 

Nondimenelonal  Transformations 

Before  developing  the  optimal  controls,  it  will  be  useful  to 
introduce  nondimens ional  variables  that  are  the  order  of  one. 

Let 


u  =  V/V  ,  x  *  gt/V 
o  ’  o 


(58) 
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The  transformed  differential  equations  are 


dx 


T  cos(a+6v> 


N  sina-siny 


(59) 


dx  “  u  {[t  8*n(a+6v)+N  cosaj coso-cosy j  (60) 

3*  =  J  CT  sin(a+8v)+N  cosal  (61) 

Consider  next  the  development  of  the  optimal  controls. 

Optimal  Controls 


The  optimal  controls  are  determined  from  the 
the  Hamiltonian  defined  by 


imlzation  of 


H  *  p  ^  +p  dl 

ru  Ht  r^/  Ht 


+P 


d* 


(62) 


where  Pu»  P^,  P^  are  costate  variables  that  must  satisfy  along  the 
optimal  trajectory 


dP 

dH 

- d  , 

dx 

du 

dP 

dH 

■  e 

dx 

'  dy 

dP. 

dH 

— 4  , 

" 

dx 

(63) 


(64) 


(65) 


Since  H  is  Independent  of  the  heading  angle  <|>,  P^  is  a  constant 

everywhere  along  the  optimal  trajectory.  Since  the  final  heading 
is  free 


(66) 


The  boundary  condition  on  the  costate  variable  corresponding  to 
the  flight  path  angle  y  is 


W 


(67) 


If  final  speed  is  free,  from  the  transversal ity  condition 


w 


(68) 


If  the  final  speed  is  prescribed 


Pu(V 


(69) 


The  optiaal  angle  of  attack  schedule  as  a  function  of  tiae  will  be 
different  for  the  free  and  specified  final  speed  cases.  Hereafter 
we  will  only  consider  the  case  where  the  final  speed  is  free.  The 
boundary  conditions  in  Equations  (66),  (67),  and  (68)  aust  be  such 
that  the  boundary  condition  in  Equation  (57)  is  satisfied.  The 
final  tiae  is  free  and  H  does  not  contain  the  tiae,  therefore 

H  *  C  (70) 

o 


everywhere  along  the  optiaal  trajectory. 

Substituting  Equations  (59)  through  (61)  and  (66)  into 
Equation  (62)  gives 


H  »  atT  -P^siny  _l 


cosy 


(71) 


where 
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(72) 


Pucos(a+6v)+  ~  P^cososin(<x+Sv) 

Am»  -P  sina+^  P  cosacosa  (73) 

N  u  u  y 

Consider  first  that  part  of  the  Hamiltonian  that  contains  the 
bank  angle 

H(o)  «  J  P^cosaCT  sln(a+6v)+N  cosa3  (74) 

There  are  only  two  choices  that  maximize  H 

a*  »  0  if  P  >  0  (75) 

"Y 

ff*  >  i  if  <  0  (76) 

The  final  flight  path  angle  is 

A 

Yltf)  »  1  if  O  =0 

(77) 

■yit^j  =  -»  if  a*  =  % 

The  optical  path  is  in  the  vertical  plane.  If  the  final  speed  is 
free,  the  optimum  is  a  half  loop  or  split-S  both  of  which  are 
illustrated  in  Figure  21.  If  the  final  speed  is  specified,  then 
an  additional  maneuver  must  be  considered.  This  third  Maneuver  is 
a  pull  up  to  stall,  a  rotation  to  a  nose  down  attitude,  and  then  a 
pull  up  to  zero  flight  path  angle  and  specified  speed.  In  the 
numerical  examples  to  be  examined  later,  we  will  consider  only 

free  final  speed.  The  only  time  that  oA  can  be  different  from  0 
or  v  is  when  there  is  an  altitude  constraint  or  the  final  heading 


is  not  «.  These  possibilities  will  not  be  addressed  in  this 
effort. 

He  can  now  drop  Equation  (61).  Equations  (72)  and  (73) 
become 


P  5 

at  *  Pucos(a+6v)+  -jjj-  sin(a+6v)  (78) 

P  6 

am  =  P„sina+  -~c—  cosa  (79) 

N  U  U 

where 

6  *  +1  for  a  pull  up  (80) 

6  -  -1  for  a  pull  down  (81) 

Consider  next  optimal  thrust  Magnitude  control. 

The  thrust  appears  linearly  in  the  Hamiltonian,  therefore, 
for  optimal  thrust  magnitude  control 


T  -  T»x 

if  XT  )  0 

(82) 

T  =  0 

if  xT  \  0 

(83) 

and  may  be  interior  if  xT*  0  for  a  finite  time.  This  situation 
will  be  examined  later  in  this  section. 

Next  consider  optimization  with  respect  to  the  thrust  vector 
angle,  6y.  Let  H(Sv>  denote  that  part  of  the  Hamiltonian  that 

contains  terms 
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H(6yJ  *  XtT 


(84) 


Haxialzatlon  with  respect  to  6y  requires  that 


d& 


«  0 


(85) 


Substitution  of  Equation  (78)  qlves 


P  S 

-Pusin(«+6v.)-f  cos(a+6v)  *  0  (86) 

The  solution  may  be  in  any  quadrant,  therefore,  the  sua  is  coa- 
puted  froa 


sin(a+6v)  *  + 


P  6 
uP 


,  COS  (01+6  )  *  + 


u 


(87) 


where 


1-&) 


Z-.  1/2 


(88) 


Substitution  of  Equations  (87)  into  Equation  (78)  gives 


*T  a  ±  P 


(89) 


Maxlaua  H(6v>  corresponds  to  aaxinua  xT,  therefore,  the  positive 
sign  holds  and 

xT  -  P  (90) 

Hv6v)  «  PT  (91) 
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■/vv’v./-.-v 


pvs  pu 

sima+6v>  =  -g-  ,  costa+6y)  *  p-  (92) 

Intermediate  thrust  can  be  ruled  out  since  this  would  require  P  to 
be  zero  which  requires  Py  and  P^  both  to  be  zero  according  to 

Equation  (88).  This  can  not  hold  for  a  trajectory  optimization 
problem,  the  costate  variables  cannot  all  be  zero.  Thus,  maximum 
thrust  magnitude  is  applied  throughout  the  trajectory. 

The  optimal  angle  of  attack  is  determined  next.  Let  H(a) 
denote  that  part  of  the  Hamiltonian  that  contains  a  terms 


Hia)  *  X^pT  +xnN 


The  necessary  condition  for  optimum  a  is 


/ 

dH\a)  d\T  dXM  dN 

da  da  A  da  "  AN  da 


0 


(93) 


(94) 


Since 


.  o 

da  d6y 


(95) 


Equation  (94)  becomes 


P  6 


P  6 


(-P.cosa-  sina)N  +(-P.  slna-t  cosa)N  *  0 
u  u  u  u  a 


(96) 


This  is  a  transcendental  equation  in  a.  Solving  for  tana  gives 
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tana  » 


(97) 


Recall  for  the  Instantaneous  turning  problem  the  optimum  angle  of 
attack  was  and  the  thrust  vector  angle  was  the  complement 
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astall  is  the  sa*e  as  the  final  condition  for  the  problem 


From  Equations  (92) 


cote  *  P, 


(102) 


A  A  A 

e  *  a  +6V 


(103) 


ibstitutlnq  this  into  Equation  (97)  results  in  a  relation  between 

and  a*  that  must  be  satisfied  everywhere  along  the  optimal 
ajectory 


cote 


-tana 


1  +  —  tana 

'■'N 


(104) 


e  relationship  between  the  optimal  angle  of  attack  and  optimal 
silon  is  presented  in  Figure  23.  Given  the  optimal  angle  of 
tack,  the  optimal  thrust  vector  angle  can  be  determined.  The 
stribution  is  presented  in  Figure  24.  At  high  angles  of  attack 
ere  CN  *  0,  Equation  (104)  reduces  to 


cote  *  -tana 


(105) 


e  solution  of  this  equation  is 


A  ir  * 
e  =  j  -a 


(106) 
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The  solution  for  optimal  6V  is  therefore 


( 107) 


This  holds  whenever  =  0.  For  small  angles  of  attack,  CN  is 

a  w 

approximately  CM  a.  Equation  (104)  is  approximately 


A  A 

e  =  2a 


(108) 


The  optimal  thrust  vector  angle  at  low  angles  of  attack  is 
approximately 


.A  A 

t>v  3  a 


(109) 


For  free  final  speed,  the  final  optimal  angle  of  attack  is  the 
stall  angle  of  attack.  The  only  way  that  higher  angles  of  attack 
can  be  reached  is  if  the  costate  function  PR  is  less  than  zero. 

This  can  be  seen  from  examination  of  Equation  (97)  or  Figure  22. 

In  addition  to  satisfying  Equation  (96),  it  is  necessary  for 
a  relative  maximum  that  H(a)  satisfy 


dzH(g) 


(110) 


Substitution  gives 


i-'H(a)  a  _x  T  N  +2  ^  ^ 

da2  T  N  ax  3<x  N  8a 


(111) 
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e  WO,  it  follows  that  X„>0.  Also,  N  >0,  N  <0,  T  >0,  N  >0 

I  N  0t~  0UX~ 

substitution  does  verify  that  Relation  (110)  is  satisfied. 

Since  the  optimal  maneuvers  are  in  the  vertical  plane,  the 
ectories  are  either  a  half-loop  or  a  split-S.  For  a  half- 
,  o  =  0 ,  6  =  1 ,  and  *  *•  For  a  split-S,  a  -  w, 

-1,  and  y(t^)  =  These  terminal  conditions  replace 

tlon  (57).  The  optimal  trajectory  is  the  maneuver  that  gives 
minimum  time  to  rotate  through  180  degrees  in  the  flight  path 
e.  This  can  only  be  determined  by  numerically  integrating  the 
oprlate  differential  equations  and  comparing  times  for  a  half- 
and  a  split-S. 

The  differential  equations  that  must  be  solved  simultaneously 
the  differential  equations  for  u,  y,  Pu,  and  P^.  Since  the 

mal  controls  are  a  function  of  PR,  the  differential  equations 

Pu  and  P^  can  be  replaced  with  a  single  differential  equation 

PR.  Differentiating  the  logarithm  of  PR  gives. 


dx 


!b  if 

u  dx 


P  dx  P  dx 
u  y 


(112) 


tituting  Equation  (59)  for  du/dx  and  expanding  Equations  (63) 
(64)  for  dPu/dx  and  dP^/dx  gives  after  some  manipulation 


R  i  i 

*  *  T  (PRcose+sine)+  —  N  ( PRsin<x-cosa) 

-  y  S(1+Pr)cosy 

tituting  Equations  (92)  for  e  gives 


(113) 
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I 


~  T  Kl +P*/1/2+  i  N  ( Pesina-cosa> -  -  6<l+P2)cosy  (114) 

U  n  U  K  U  i\ 


The  right  side  is  a  function  of  the  costate  variable  PR,  the  state 

variables  u  and  -y,  and  the  control  variables  T  and  «.  The  dif¬ 
ferential  equations  for  u  and  y  are  a  function  of  the  state  and 
control  variables.  The  control  variable  a  is  a  function  of  the 
costate  variable  and  the  thrust  is  maximum.  Thus,  three  differen¬ 
tial  equations  result,  one  each  for  PR,  u,  and  y.  The  only 

unknown  is  the  initial  value  for  PR.  But  the  final  value  is  known 

from  Equation  (68).  Thus,  a  single  parameter  search  on  the  ini¬ 
tial  value  of  PR  is  required.  An  initial  estimate  is  made,  the 

three  differential  equations  are  integrated  forward  in  time,  and 
the  optimal  trajectory  is  that  trajectory  that  gives  =  0. 

Later  on  in  this  section,  we  will  develop  some  numerical  results. 
Consider  next  the  problem  where  the  thrust  vector  angle  is  fixed 
at  zero. 

Zero  Thrust  Vector  Angle 

Equations  (92)  do  not  apply.  Equation  (94)  for  optimal  a 
becomes  upon  substitution  of  Equations  (78)  and  (79)  with  Sv*  0 


dH(g) 

da 


V  -atn 


+A 


N 


£N_ 

da 


0 


(115) 


This  equation  is  different  from  Equation  (94)  because  of  the  x^T 
term.  Equation  (110)  becomes 


d?H(ou 

aa' 


-\  T  N  -2X  ~  +X  —  — 
\TL  *nn  zxt  aa  +xN  2 

a(X 


0 


(116) 
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Solving  Equation  (115>  for  xw  and  substituting  gives 


.  2 

oa 


-A.n 


T  +2 


^  +|n  -  *-*  | 

da  1  «>a2  ] 


iili 


w 


T  +r^ 
oa 


v  117 ) 


Since  the  first  derivative  of  N  is  non  negative  and  the  second 
derivative  is  non  positive ,  the  bracket  sum  is  non  negative  and  it 
follows  that  if  Relation  (117)  is  satisfied,  then 

at  ,  0  (118) 

The  thrust  magnitude  is  maximum  except  for  possible  Interior 
thrust.  Consider  this  possibility.  For  interior  thrust,  xT  *  0 

for  a  finite  tine,  thus  Equation  (115)  reduces  to 

■  *n(t  +  £)  *  0  <“9» 

The  only  possibility  is 

*N  »  0  *120) 

If  both  aw  and  xT  are  zero,  then  Equations  (78)  and  (79)  reduce  to 
P  6 

P^cosa*  — sina  *  0  (121) 

P  6 

cosa-P^sina  =  0  (122) 


These  two  equations  combined  together  give 


*  *  fr) 


*  u 


t  123) 


This  can  hold  if  and  only  if 


P  =  P  *  0 
u  y 


*  124 ) 


This  can  not  hold  for  the  optimal  solution.  Therefore  xT#  0  for  a 

finite  time  and  maximum  thrust  is  applied  everywhere  along  the 
optimal  trajectory. 

The  optimal  angle  of  attack  is  the  solution  of  Equation 
(115).  Substituting  Equations  (78)  and  (79)  gives 


P  6 


P  6 


<-Psina+  — cosaMT  +N  ) -(P  cob<x+  — *-  sina)N  *  0  (125) 

u  u  a  u  u 


Since  the  optimal  angle  of  attack  is  between  zero  and  ninety 
degrees.  Equation  *125)  can  be  rewritten  as 


tana  * 


T  j  Na  -  PRN 

Pc  *T  +  N  )  +  N 

(X 


*126) 


The  optimal  angle  of  attack  is  the  solution  of  Equation  (126). 


At  the  final  time,  Pu*  0  from  Equation  *64).  Substitution 


into  Equation  (126)  gives 


tana 


T  +N_ 


W 


*127) 


Comparison  with  Equation  *26)  shows  that  the  final  value  of  the 
optimal  angle  of  attack  is  equal  to  the  optimal  instantaneous 
value.  Thus,  if  the  final  speed  is  higher  than  the  critical 


33 


speed,  the  final  angle  of  attack  is  close  to  <*stall*  The  only 

time  that  high  angle  o£  attack  could  occur  is  at  speeds  less  than 
the  critical  speed,  in  air  combat,  this  situation  is  to  be 
avoided  since  the  engangement  is  similar  to  that  between  helicop¬ 
ters  and  fighters.  But  from  strictly  a  single  vehicle  minimum 
time  turn  problem,  the  angle  of  attack  could  attain  large  values 
if  the  speed  is  very  low. 

Numerical  Examples 

Optimal  controls  tor  zero  thrust  vectoring  will  not  be  em¬ 
phasized  since  it  was  shown  in  the  last  section  that  high  angle  of 
attack  maneuvers  can  result.  With  thrust  vectoring  there  is 
uncertainty  as  to  whether  or  not  PSM  is  optimal.  At  this  point 
all  that  we  know  is  that  the  final  value  of  the  optimal  angle  of 
attack  equals  the  stall  value. 

Optimal  and  nominal  trajectories  were  developed  for  Initial 
conditions  corresponding  to  0.6  Mach  Number  and  20,000  feet 
altitude.  The  nominal  trajectories  correspond  to  constant  angle 
of  attack  and  thrust  vector  angle.  The  angle  of  attack  was 
“stall '  Optimal  trajectories  were  in  the  vertical  plane.  The 

pull  up  maneuver  required  less  time  to  rotate  through  180  degrees 
than  did  a  split-S  maneuver. 


The  numerical  results  are  presented  in  Table  1.  Six  dif¬ 
ferent  cases  are  presented.  The  thrust  to  weight  ratio  of  0.6  is 
a  typical  value  for  the  F15  at  20,000  feet  and  0.6  Mach  Number. 
Higher  values  of  T  were  introduced  to  determine  the  effect  of 
increasing  thrust.  The  first  two  cases  were  nominal  trajectories 


flown  at  C, 


i .  e , 


max 


“stall*  dif^erence  was  *n  the  thrust 


vector  angle.  Thrust  vectoring  reduced  the  time  by  approximately 
three  and  one  half  seconds.  The  last  four  cases  are  optimal 
results.  The  third  and  fourth  cases  show  the  difference  between  a 
[ull  up  and  a  split-S  maneuver.  The  pull  up  is  better. 
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Comparison  of  the  third,  fifth,  and  last  case  shows  that  the  tine 

decreases  as  the  thrust  to  weight  increases.  This  is  consistent 

with  the  results  from  the  previous  section.  There  the  turning 

rate  increased  with  increasing  thrust  to  weight  ratio.  The  angle 

of  attack  did  reach  values  where  Cu-  C.,  .  The  thrust  vector 

N  Nnax 

angle  went  to  ninety  degrees. 


TABLE  1 

TIME  TO  TURN  RESULTS 


T/W 

MANEUVER 

MAX  a  (Deg) 

MAX  6v  (Deg) 

TIME  (SEC) 

0.6 

PULL  UP 

astall 

90-a  .  , . 
stall 

10.8 

0.6 

PULL  UP 

“stall 

0 

14.5 

0.6 

PULL  UP 

50 

90 

9.6 

0.6 

SPLIT-S 

30 

60 

12.3 

1.0 

PULL  UP 

55 

90 

7.9 

1.2 

PULL  UP 

60 

90 

7.3 

Based  upon  the  limited  number  of  numerical  results  examined 
here,  the  conclusion  is  that  for  minimum  time  turns,  post  stall 
maneuvering  is  optimal.  This  holds  for  both  thrust  vectoring  and 
nonvectoring. 

In  Figures  25  through  29,  u,  y,  PD,  ot,  and  S„  are  presented 
as  a  function  of  the  nondimensional  time  and  for  initial  values  of 
Pg  equal  to  0.3,  0.4,  0.5,  0.6,  and  0.7.  These  results  correspond 
to  the  third  case  in  Table  1.  The  trajectories  for  Pg  (0)  =  0.3 
and  0.4  do  not  generate  enough  lift  to  fly  through  the  vertical, 
thus  the  heading  can  not  change  by  180  degrees.  The  optimal 
trajectory  corresponds  approximately  to  PR  (0)  =  0.5.  In  Figure 
28  it  can  be  seen  that  the  maximum  angle  of  attack  is  ap¬ 
proximately  50  degrees.  From  Figure  25  it  can  be  observed  that 
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the  final  speed  is  very  low.  These  results  demonstrate  that  high 
angle  of  attack  maneuvers  are  optimal  with  thrust  vectoring. 

Sumfary 

Minimum  time  turns  in  the  subsonic  region  require  maximum 
thrust  throughout  the  maneuver.  This  holds  for  both  vectored  and 
nonvectored  thrust.  For  the  vectored  thrust  case  the  final  angle 
of  attack  equals  astail  which  matches  the  instantaneous  problem 
which  holds  everywhere.  For  nonvectored  thrust  the  final  angle  of 
attack  equals  that  for  maximum  instantaneous  turning  rate. 

Optimal  maneuvers  employ  PSM.  As  the  thrust  to  weight  ratio 
Increases  so  does  the  maximum  angle  of  attack  along  the 
trajectory. 
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SECTION  V 


CONCLUSIONS 

/ 

'  ^'The  purpose  of  this  effort  was  to  determine  whether  or  not 
there  are  any  apparent  performance  Improvements  through  Post  Stall 
Maneuvers  iPSM)  or  thrust  vectoring.  PSM  is  a  result  of  high 
angles  of  attack,  greater  than  the  stall  value.  Two  different 
problems  were  addressed.  The  first  examined  instantaneous  turning 
performance.  The  second  focused  on  minimum  time  turn  problems. 

The  impact  of  both  vectored  and  nonvectored  thrust  was  considered. 
It  was  proven  that  minimum  time  turns  fall  in  the  vertical  plane. 

For  maximum  instantaneous  turning  rate  and  nonvectored  thrust 
high  angle  of  attack  results  if  the  speed  is  less  than  the  criti¬ 
cal  speed  which  is  approximately  Mach  0.2.  The  optimal  angle  of 
attack  approaches  ninety  degrees  as  the  speed  approaches  zero. 

For  vectored  thrust,  the  stall  angle  of  attack  is  optimum  and  the 
thrust  vector  angle  i3  the  complement  of  the  angle  of  attack.  PSM 
is  not  optimal  if  thrust  vectoring  is  available. 

For  minimum  time  turning  performance,  maximum  thrust  is 
optimal.  At  the  end  of  the  trajectory,  the  optimal  angle  of 
attack  equals  that  tor  maximum  instantaneous  performance.  PSM  is 
optimal.  As  the  thrust  to  weight  ratio  increases,  the  maximum 
angle  of  attack  increases,  j. 

Future  studies  should  address  the  following  aspects: 

1.  Unsteady  aerodynamics 

2.  Analysis  of  aircraft  pointing 

3.  Compare  turning  versus  pointing 

4.  Air  combat  analysis  MvN 
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FIGURE  4.  Steady  State  and  Dynamic  Lift 
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Altitude  =  20,000  feet 


FIGURE  6.  Turning  Rate  as  a  Function  of  Angle  of  Attack  and  Mach  Number  Variations 
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FIGURE  7.  Turning  Rate  as  a  Function  of  Angle  of  Attack  and  Mach  Number 


f  Attack  as  a  Function  of  Controls  and  Mach  Number 


FIGURE  9.  Turning  Rate  as  a  Function  of  Controls,  Altitude,  and  Mach  Number 


Altitude  =  20,000  feet 


Function  of  Thrust  Variations 


Altitude  =  20,000  feet 


FIGURE  11.  Optimum  Angle  of  Attack  as  a  Function  of  Thrust  Variations 


Altitude  =  20,000  feet 
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FIGURE  12.  Maximum  Turning  Rate  as  a  Function  of  Aerodynamic  Variations 


Altitude  =  20,000  feet 


FIGURE  13.  Optimum  Angle  of  Attack  as  a  Function  of  Aerodynamic  Variations 
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FIGURE  14.  Maximum  Turning  Rate  as  a  Function  of  Thrust  Vector  Angle 


FIGURE  15.  Angle  of  Attack  for  Maximum  Turning  Rate  as  a  Function  of  Thrust  Vecto  Angle 


Altitude  =  20,000  feet 


FIGURE  16.  Sustained  Turning  Rate  as  a  Function  of  Thrust  Vector  Angle 


FIGURE  17.  Angle  of  Attack  for  Sustained  Turning  Rate  as  a  Function  of  Thrust  Vector  Angle 
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FIGURE  18.  Maximum  Turning  Rates 


FIGURE  19.  Optimal  Angles  of  Attack 
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FIGURE  22.  Optimal  Angle  of  Attack  as  a  Function  of  the  Costate  Function 


FIGURE  23.  Optimal  Angle  of  Attack  as  a  Function  of  Optimal  Epsilon 
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FIGURE  26.  Flight  Path  Angle 


20,000  feet 


0.0  0.1  0.2  0.3  0. 

NONDIMENSIONRL  TIME 


1*1* 


28. 


to 

o 


67 


U.S.Qo»«rnm«nt  Printing  Office  1985  —  559-065/20856 


•  v  /-  .t  .-  V- 


*•  •  »  , '«  *•  . ’►  .  \  *»  *,  \  * ,  8  \  •  8.  '  »  *  8  *  «  '  *  ■* 


j 


FIGURE  29.  Thrust  Vector  Angle 


